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We investigate the entanglement spectrum near criticality in finite quantum spin chains. Using 
finite size scaling we show that when approaching a quantum phase transition, the Schmidt gap, 
i.e., the difference between the two largest eigenvalues of the reduced density matrix Ai, A2, signals 
the critical point and scales with universal critical exponents related to the relevant operators of the 
corresponding perturbed conformal field theory describing the critical point. Such scaling behavior 
allows us to identify explicitly the Schmidt gap as a local order parameter. 
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The characterization of phase transitions in strongly 
correlated systems has been traditionally based on the 
behavior of expectation values of local operators and two 
point correlators. This is the spirit of the "standard" 
Ginzburg-Landau scenario of second order phase transi- 
tions where many-body systems order (locally) for given 
values of a "control" parameter g of the Hamiltonian. 
The order, which is associated to the breaking of some 
symmetry, is manifested by an order parameter, Q, whose 
expectation value in the ordered phase is different from 
zero while it becomes exactly zero at criticality, i.e., for 
g = gc- Near criticality, order parameters and other phys- 
ical quantities describing the many-body system exhibit 
scaling behavior, i.e., a power dependence with \g — gd- 
The exponents of these power laws are called critical ex- 
ponents. The scaling exponents reflect the universality 
class of the theory, which is independent of the micro- 
scopic details of the model but depends only on global 
properties such as the symmetries and dimensionality of 
the Hamiltonian. 

Recently, an alternative approach to understanding 
quantum many-body systems and their simulatability [H , 
exploits the entanglement content of the ground states 
of such systems. First attempts within this approach fo- 
cused on the study of quantum phase transitions (QPT) 
in spin chains 0, 0] ■ A good figure of merit to measure 
the bipartite correlations embedded in the ground state 
of a spin chain of length L, is the entanglement entropy 
(EE). It refers to the von Neumann entropy of a block 
of £ < L spins, S ~ — ^^A^logAi, where the Ai, are 
the eigenvalues of the reduced density matrices obtained 
after the bipartite splitting of the ground state. Its fi- 
nite size scaling, i.e., the dependence of the EE with the 
size of the block i, shows remarkable properties At 
criticality, the EE diverges logarithmically as S" ~ c log £ 
where c is the so-called central charge of the cor- 
responding QPT as provided by conformal field theory 
(CFT). This universal logarithmic behaviour of the en- 



tanglement entropy at criticality underpins the conformal 
invariance of QPT in ID, and leads also to a universal 
-depending asymptotically only on c- distribution of the 
eigenvalues of the reduced density matrix [9|. On the 
other hand, outside but close to criticality, when the sys- 
tem is gapped, the EE scales as S" ~ clog(^) being 
^ the correlation length that sets the relevant scale for 
long-distances physics. It is in this regime when the sys- 
tem is not any longer entirely constricted by the central 
charge, and its description depends also on the full set 
of Hamiltonian parameters tuned to move away from the 
critical point. 

Thus, it is to be expected that away from criticality, 
further information about the many-body system that is 
not included in EE can be obtained from the set of pa- 
rameters {Ai} 13, Hi- We shall denote this set as the 
entanglement spectrum. Indeed, it is known that the 
(topological) Haldane phases appearing for integer spin 
chains arc characterized by a double (or higher) degener- 
acy of the entire entanglement spectrum, resulting from 
applying the symmetries of the Hamiltonian to the eigen- 
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states of the reduced density matrix 

Here, we examine the behavior of the entanglement 
spectrum for finite spin chains, of length L, in the vicin- 
ity of a QPT. For simplicity, we assume a symmetric 
bipartite splitting and we show that when approaching a 
QPT, the difference between the two largest (nontrivially 
degenerate) eigenvalues of the reduced density matrix 
Ai, A2 signals correctly the critical point and scales with 
the universal critical exponents. We relate this result to 
the relevant operators of the corresponding CFT valid to 
describe perturbations from the critical QPT point. 

Before proceeding further towards our main result, let 
us just clarify the notation used throughout. For any 
bipartite splitting, the ground state of the spin chain 
can be expressed by its Schmidt decomposition \ipGs) = 
J2i I'^f )® I'/'f ) ) where Ai > are the Schmidt eigen- 
values with respect to the partition (left or right) sorted 
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in decreasing order and \<t>f) and |0f ) are the Schmidt 
eigenvectors [l^. Our study refers to partitions of the 
spin chain in real space and it differs from previous stud- 
ies in the quantum HaU regime and some spin chains 
lol 14 1 , in which the entanglement spectrum is analyzed 
after a cut in momentum space is performed. 

We first review the finite size effects exhibited by the 
entanglement spectrum at criticality. In the framework 
of CFT (see [3 US] and references therein) for a chain 
of size L with periodic boundary conditions, the reduced 
density matrix of a block of size £ can be expressed as 
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where c is the central charge, Lq is the zero-generator 
of the chiral Virasoro algebra, Zg{q) = Tr[g^''] = 1 + 
nig"i + ^2?"^ + ... is the partition function of the sub- 
system where q = exp(i27rr) with r = {iK)/[\og{i/r])], -q 
is a regularization cutoff and k a positive constant. The 
Schmidt gap, defined as AA = Ai — A2, can be expressed 
as 
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where, in the last term, we made use of the fact that, at 
criticality, the largest eigenvalue of the reduced density 
matrix, the so-called single copy entanglement , corre- 
sponds to half of the von Neumann entropy Ai ~ t~'^l^'^ . 
The coefficients on > 0, with degeneray n^, related to the 
scaling dimension of the relevant operators of the theory, 
correspond to the eigenvalues of the operator Lq. There- 
fore, from CFT arguments, for finite systems at criticality 
the Schmidt gap is solely determined by c and the small- 
est eigenvalue ai. Notice that in the thermodynamic 
limit the Schmidt gap closes. 

Here we study the entanglement spectrum {A^} in the 
vicinity of a QPT by varying a single parameter of the 
Hamiltonian describing the spin chain. In many known 
cases, this ensures an unambiguous identification of this 
operator with the operator content of the CFT. In these 
cases, the mass gap and corresponding order parameter 
scale as a power-law in |g — 5c| — ^ 0, with certain critical 
exponents. The latter ones depend on the central charge 
and on the primary operator perturbing the CFT. 

Our main result shows that the Schmidt gap develops 
the same scaling behavior as the one found for the mass 
gap and the magnetization, when approaching a QPT. 
That is, when properly scaled, the quantities AA(L, g) 
obtained for different chain lengths L can be made to 
cross in a single point — which signals gc — and to col- 
lapse into a single curve, with the same universal scaling 
exponents of the QPT. It is indeed remarkable that the 
information encoded in just these two eigenvalues repro- 
duces global properties of the ground state out of criti- 
cality, where equation (2) is not valid anymore. 
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FIG. 1: Transverse Ising model. Top-Left: Schmidt gap at the 
critical point ( J = B^) as a function of the chain size. The nu- 
merical data (symbols) are compared to a fitting from Eq. ([2]) 
(dashed) and to a simple power law (dotted). Top-Right: En- 
tanglement spectrum: the first 8 eigenvalues Ai are shown as 
a function of J/Bz for L = 12000. Bottom-Left: Schmidt gap 
as a function of J/B^ for different lengths. Bottom-Right: 
ESS analysis of AA for L = 768, 1536, 3072, 6000, 12000. 



We illustrate our results by analyzing first the entan- 
glement spectrum and the Schmidt gap in the integrable 
transverse field Ising model. We move then to the longi- 
tudinal spin-1/2 Ising model and the nonintegrable uni- 
axial spin-1 chain which we solve numerically. Let us 
stress that at criticality, and for sufficiently long chains 
our results agree with those obtained by CFT. 

Spinl/2 Ising model- Wc start with the following 
spin-1/2 Ising model: 
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where ak=x,y,z denotes the Pauli matrices. We assume 
open boundary conditions. For = 0, Eq. ^ reduces 
to the paradigmatic transverse-field Ising model with a 
critical point at J = Bz and central charge c = 1/2. 
To compute the entanglement spectrum we employ the 
Jordan- Wigner and Bogoliubov transformations to map 
the transverse Ising model into a system of noninteract- 
ing fermions. We then proceed as in [l^ to obtain the 
Schmidt eigenvalues Ai of the reduced density matrix af- 
ter cutting a chain of size L (from 48 to 12000) in two 
symmetric halves {£ = L/2). Our results at criticality, 
J = Bz, are summarized in the top- left panel of Fig. [U 
where AA, shown as a function of L, is fitted by a simple 
power law function L~'^ and by the expression given by 
Eq. ([2]) using c = 1/2, ai = 1/8, and 77 and k as fitting 
parameters. Both fittings are very accurate but when 
looking at the marginals of the fitting functions, Eq. 
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matches better the Schmidt gap values, in particular for 
the largest chain size, i.e. L = 12000. 

For J > Bz, the order parameter = L^^J2iWx)^ 
scales, in the vicinity of the critical point, as Mr ~ 
\J/Bz — with critical exponent /3 = 1/8 while the 
correlation length, ^ ^ \J/Bz — l]""^ scales with v ~ 1. 
We study now perturbations from the critical point by 
changing the parameter J/Bz while keeping Bx = 0. 
The results for the entanglement spectrum are plotted 
in Fig. [T] (Top- Right panel), showing a tendency of the 
eigenvalues to collapse at the critical point (notice the 
vertical logarithmic scale). The doublets appearing for 
J > Bz are due to the unbroken Z2 symmetry of the 
ground state. The behaviour of the Schmidt gap AA 
as a function of J/Bz for different lengths is shown in 
FigH] (Bottom-Left panel). The Schmidt gap closes very 
rapidly when J approaches Bz- For finite size systems, 
the asymptotic behavior corresponding to the thermo- 
dynamic limit can be retrieved using finite size scaling 
(FSS) ■ Near criticality, the dependence of the order 
parameter Q with the finite size L is given by 
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where v characterizes the divergence of the correlation 
length while /3q is the order parameter critical expo- 
nent: Q ~ I5 — ffd^"- To obtain critical exponents, we 
change the coefficients /ii_2 until we observe the collapse 
of AA L'^i as a function of \g — gc for all used lengths 
L in the simulations. In this way we extract /3q = /ii//i2 
and V = 1///2- The results of the FSS of the Schmidt gap 
are depicted in Fig. [1] (Bottom- Right panel). From the 
scaling of the Schmidt gap we obtain the critical expo- 
nents /3aa = 0.124± 0.002 and i^aa = l.OOiO.Ol. These 
values coincide with the critical exponents of the Ising 
universality class, suggesting that the Schmidt gap in the 
vicinity of the critical points scales universally with the 
same critical exponents. We remark that FiglT] contains 
results computed using the symmetric superposition of 
the two degenerate ground states in the ferromagnetic 
phase as obtained through the Jordan- Wigner transfor- 
mation. The same results for the critical exponents arc 
obtained if the Z2 symmetry is explicitly broken and only 
one of the two degenerate ground states is considered. 

The surprising, yet interesting, results obtained for the 
transverse Ising model pushes us to look at the same 
critical point of the transverse Ising model J = Bz but 
adding also a very weak longitudinal field perturbation, 
i.e. B^ < Bz in Eq. This integrable model [ISl 

has the same critical point as the transverse Ising model, 
Bx = 0, and thus a central charge c = 1/2, but belongs 
to a different universality class. This implies for instance 
that its critical exponents {v = 8/15 and /3 = 1/15) do 
not coincide with the exponents of the transverse Ising 
Hamiltonian. However notice that the ratio /3/i^ = 1/8 is 
the same in both cases as predicted by scaling dimen- 
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FIG. 2: Ising model with longitudinal perturbation, DMRG 
calculations. Left: The first 6 Schmidt eigenvalues as a func- 
tion of B^l in logarithmic scale for L = 384. The critical 
point corresponds to Ex — >■ 0. Right: FSS analysis of AA for 



sional analysis. We use the density matrix renormal- 
ization group (DMRG) j2l| to find the ground state of 
the Hamiltonian Eq. ^ in the vicinity of the critical 
point and to obtain its entanglement spectrum. The re- 
sults are shown in Fig. [2l In the left panel we show the 
first 6 eigenvalues \i which, when approaching the crit- 
ical point, tend to form a band as expected from GET 
22I |23|. In the right panel we show the FSS analysis 
of the Schmidt gap AA from which wc infer the esti- 
mate for the critical exponents: v = 0.50 ± 0.05 and 
/?AA = 0.055 ± 0.005 which arc in good agreement with 
the exact values of the corresponding exponents, con- 
firming that the scaling exponents of the Schmidt gap 
properly reflect the perturbation on the critical point. 

Finally, we consider also the nonintegrable spin-1 
Heisenberg chain with uniaxial anisotropy which presents 
a very rich phase structure The Hamiltonian de- 

scribing such a system is given by 
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where Si = {Sxi, Syi, Szi) are the zth site angular momen- 
tum operators and Hi[9) = cos{6)Si ■ Si+i + sin{9){Si ■ 
Si+i)'^. For D = 0, the model is the well known 
bilinear-biquadratic spin-1 chain, whose phase diagram 
as a function of 9 G [— tt; tt] is well established. For 
— 7r/4 < 9 < 7r/4, the system is in the aforementioned 
Haldanc phase. A sketch of the phases surrounding the 
Haldane phase when the uniaxial anisotropy D is present 
is depicted in the inset of Fig. [3l We concentrate our 
analysis here in the Neel-Haldane phase transition along 
9^0 using the uniaxial anisotropy, D, as the control 
parameter but the same results are obtained if the tran- 
sition is crossed at a fixed value of D and using 9 instead 
as the control parameter. In the Neel phase, the stag- 
gered magnetization per site Mz = J2i{~^y i'^zi) is 
the order parameter. As before, we calculate the ground 
state and the Schmidt eigenvalues after a bipartite cut 
of the spin chain into two symmetric halves, obtain the 
Schmidt gap and perform FSS for different chain sizes. 



4 




' ■ ■ ■ ' u. / 

-0.33 -0.31 -0.29 -0.27 -0.25 -2 -1.5 -1 -0.5 0.5 1 1.5 2 

D (-0.315- D)*L°-^'' 



FIG. 3: Spin-1 model. Left: The entanglement spectrum 
for 6 = and as a function of D nearby the Neel-Haldane 
quantum phase transition (for clarity only the first 6 eigen- 
values are depicted). Dotted, dashed and continuous lines 
correspond to spin chains of length L=192,384 and 768, re- 
spectively. The location of the critical point obtained after 
FSS analysis of the Schmidt gap is D = -0.315. Right: FSS 
analysis of the Schmidt gap together with a schematic plot of 
the critical points marking the transition between Neel (N), 
Haldane(H) and Large-D (L) phases for 6* = as a function 
of D. 



Observable 


Dc 


/? 


V 


M, (QMC) 


-0.316 


0.147 


1.01 


M, (DMRG) 


-0.315 


0.11 


1.01 


AA (DMRG) 


-0.315 


0.11 


1.04 



TABLE I: Neel-Haldane transition in the spin-1 model. The 
first line contains results obtained with quantum Monte Carlo 
of Ref. [2^. The other two lines contain the FSS anal- 
ysis from the order parameter (staggered magnetization) 
and the Schmidt gap respectively for system sizes of L = 
96, 192, 384, 768. 



In our DMRG simulations, we apply a small magnetic 
field to the first spin to select one of the two degenerate 
ground states. This technical trick has the advantage of 
stabilising and speeding-up the algorithm. Notice that 
in the thermodynamic limit only one of the two degen- 
erate ground states is selected: our procedure allows to 
investigate correctly this limit. The entanglement spec- 
trum is shown in the left-panel of Fig. [3] where the ex- 
pected doubly degeneracy in the Haldane phase is evi- 
dent. The FSS analysis of the Schmidt gap is shown in 
Fig. [3] (right-panel). Our estimates for the location of 
the critical point Dc and critical exponents /?, v obtained 
using FSS on the staggered magnetization and on the 
Schmidt gap are summarized in Table H] and compared 
with those obtained using quantum Monte Carlo calcu- 
lations [2^. The critical exponents obtained from the 
Schmidt gap confirm the well known fact that the Neel- 
Haldane transition belongs to the same universality class 
as the transverse Ising model. There is a 12% discrep- 
ancy in the critical exponent ji that can be attributed to 
the relatively small size of the chains, as compared to the 
12000 sites used in the transverse spin-1/2 Ising model. 
Summarizing, we have analyzed the entanglement 



spectrum {A^} for finite size systems near the QPT of 
some integrable and nonintegrable spin models. Notably, 
the information encoded just in the difference between 
its two largest eigenvalues (Schmidt gap), reproduces the 
scaling behavior of conventional order parameters closing 
at criticality and displaying universal critical exponents 
on its finite size scaling. The first result stems naturally 
from the fact that at criticality and in the thermody- 
namical limit the entanglement spectrum must become a 
continuous distribution of eigenvalues implying the clo- 
sure of the Schmidt gap. While finite size effects generally 
prevent such closure, the presence of a symmetry which 
is preserved in the Schmidt decomposition of the ground 
state, like the Zi parity in the transverse Ising or in the 
considered transitions in the spin-1 chain forces the 
entanglement spectrum to organize in multiplets making 
finite size effects much less important. The scaling be- 
havior of the Schmidt gap when approaching a QPT is 
instead less trivial and points out that different models 
belonging to the same universality class when perturbed 
with a single operator share a common Schmidt scaling 
behavior, while fulfilling majorization relations [lij . Al- 
though we have not found such scaling for other differ- 
ences between Schmidt eigenvalues, we cannot rule out 
the finite size effects and leave this question for a forth- 
coming study [2^. 
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